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Abstract Multiple rebrightenings have been observed in the multiband afterglow of GRB 
030329. Especially, a marked and quick rebrightening occurred at about t ^ 1.2 x 10^ 
s. Energy injection from late and slow shells seems to be the best interpretation for these 
rebrightenings. Usually it is assumed that the energy is injected into the whole external 
shock. However, in the case of GRB 030329, the rebrightenings are so quick that the 
usual consideration fails to give a satisfactory fit to the observed light curves. Actually, 
since these late/slow shells coast freely in the wake of the external shock, they should 
be cold and may not expand laterally. The energy injection then should only occur at 
the central region of the external shock. Considering this effect, we numerically re-fit the 
quick rebrightenings observed in GRB 030329. By doing this, we were able to derive 
the beaming angle of the energy injection process. Our result, with a relative residual of 
only 5% — 10% during the major rebrightening, is better than any previous modeling. 
The derived energy injection angle is about 0.035. We assume that these late shells are 
ejected by the central engine via the same mechanism as those early shells that produce the 
prompt gamma-ray burst. The main difference is that their velocities are much slower, so 
that they catch up with the external shock very lately and manifest as the observed quick 
rebrightenings. If this were true, then the derived energy injection angle can give a good 
measure of the beaming angle of the prompt 7-ray emission. Our study may hopefully 
provide a novel method to measure the beaming angle of gamma-ray bursts. 
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1 INTRODUCTION 

GRB 030329, wh ich is one of the brightest gamma-ray burs ts (GRBs), has a fluence of about 1.18 x 
10^^ ergs- cm~^ dRicker et al.ll2003t I Vanderspek et al.ll2004l) . It is also very close to us. The redshift is 



z = 0.168 dGreiner et alj 2003|). At the same tim e, GRB 030329 is u nambiguously c onfirmed to be as 



sociated with supernova ( Hiorth et al.l l2b03; Stan ek et"ail l2003; Math eson et al.ll2003l) . Because of these 



important characteristics, this GRB has attracted a lot of attentions. Numerous and detailed multiband 
afterglow observations have been accumulated. Multiple rebrightenings were observed in the afterglow 
of GRB 030329. Especially, a quick and marked rebrightening occurred at about t ^ 1.2 x 10^ s. 

Several models have been proposed to explain the quick and marked rebrightenings. (i) Density- 
jump model: if the external shock encounters a sudden densi ty variance of the ci r cum-burst medium , 
the emission of the afterglow may be enhanced temporarily jWang & Loebll2000t iLazzati et al.ll2002l ; 
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lNakar&Pirarill2003t iDai & Wul I200I iTam et all l2005h . But for GRB 03032 9, it has be en argued 
that the density-jump model could not produce the obse rved rebrightenings jNakar & Pi ran 2003] 
IWilUngale et"aLll2004 

iPiran et all 120041; iHuang et al.ll2006|). (ii) The two-component jet model is an- 
other kind of choice ( Berger et al. 120031) . But Huang et al. ( 20061) made a detailed numerical calculation 
and found that although the two-component jet model could bas ically reproduce t he ov erall R-band 
light curve, it was unable to explain the steep rebrightenings. (iii) IWillingale et al] ( |2004|) discussed a 
possibility that the rebrightening might be due to a supernova component, but their model is unlikely to 
produce multiple rebrightenings. 

Actually, the energy-injection model seems to be the best explanation for the steep rebrighten- 
ings ([Huang et al 2006). Granot et al. (2003) have made a general analysis on this interpretation and 
iHuang et^] ( f2006l) presented a detailed numerical study. In Huang et al.'s study, it is interesting to note 
that although the fitting result was much better than previous models, the theoretical rebrightenings still 
could not be as steep as observations. This problem was probably caused by their simple assumption 
that the energy carried by the late (and initially slower) shells was injected into the whole external shock 
homogeneously. In reality, it was probable that the late/slow shells, which coast freely in the wake of 
the pre vious external shock, should be cold and do not expand laterally, as illustrated bv lOranot et al.l 
(l2003h . So when the late/slow shells catch up with the external shock, the energy should only be injected 
into the central region of the external shock. We call this scenario as a localized energy-injection sce- 
nario. We conjecture that the energy-injection angle can be derived from the rapidness of the observed 
rebrightenings. Since these late/slow shells basically may have the same origin as those early shells that 
produce the internal shocks and the main GRB, we argued that this method could be used to hint the 
degree of collimation in GRBs. This will be a novel way to measure the beaming angle of GRBs, in 
addition to the traditional jet-break timing method. 

In this paper, we will try to reproduce the multiband afterglow light curves of GRB 030329 nu- 
merically, based on the localized energy-injection scenario. The structure of our paper is organized as 
follows. Section 2 is a detailed description of the model. In Section 3, we study the effects of various 
parameters on the rapidness of the rebrightening. In Section 4, we simulate the observed multiband after- 
glow light curves of GRB 030329 in the framework of our model. Section 5 presents a brief discussion 
and the conclusion. 



2 LOCALIZED ENERGY-INJECTION 

According to the standard fireball model, when the fireball plows throu gh the circum-burst medium, it 



will produce a strong blastwave that accelerates the swept-up electrons ( Piran|[l999l : Ivan Paradijs et al.l 
HOOO; Meszaros 2002; Zhang & Meszaros 2004; Piran 2005). The afterglow is produced by these accel- 
erated electrons. In this study, we will only consider the synchrotr on radiation frorn electrons, althoug h 
inverse Compton scattering may also play a role in some cases dWei & LulbOOOl: ISari & Esinll200lh . 
The detailed condition of afterglow is quite complicated. For instance, the blastwave may evolve from 
highly radiative regime (the radiation efficiency e = 1) to adiabatic regime (e = 0). The blastwave will 
experience the transition from ultrarelativistic phase to Newtonian phase. Meanwhile, the effect of lat- 
eral expansion of the jet, the conditi on of the circum-burst medium (homogeneous or wind- like) and 
the equal arrival time surface effect (IWax marilll997t ISarilll998t iPanaitescu & Meszaros|[l998l) are also 
important factors which influence the evolution of the afterglow. In order to reproduce the complicated 
behavior of the afterglow of GRB 030329, all these factors s hould be considered . For this purpose, nu- 
merical evaluation is more efficient than analytical method. [Huang et"ail d 19991) and iHuang & Cheng[ 
([2003) have developed an effectiv e method which can ea sily take into account all the complex con- 
ditions addressed above. Recently. Ivan Eerten et al. I ( I2OIOI) developed a more accurate code. But since 
Huang et al.'s code is more convenient for numerical solution, we will use their code for the current 
study. 

To fit the over all multiband afterglow of GRB 030329, we basically need a two-component jet 
(iBerger et al.ll2003[) : a narrow jet with a small opening angle (0o,n) and a very large initial Lorentz 
factor (70,11), and a wide jet with a relatively larger opening angle (^^o,w) but with a smaller initial 
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Lorentz factor (70, w)- These two components are coaxial, and both of them may expand laterally when 
they push through the circum-burst medium. The narrow jet can produce the observed prompt GRB 
emission and the early afterglow, while the wide jet can contribute significantly to the late afterglow. 
In our localized energy-injection model, the key ingredient is the late/slow shells that play the role of 
energy-injection. We conjecture that these shells should be launched by the central engine via the same 
mechanism as the narrow jet component. The only difference is that they are launched slightly later, and 
most importantly, with much lower velocities. They catch up with the external shocks very lately, so 
that they could only manifest as energy-injections and rebrightenings in the afterglow, but not prompt 
gamma-ray emission. It is then reasonable to assume that the half opening angle of these late/slow shells 
should be similar to that of the narrow jet, i.e. (?o,n- 

Before encountering the external shock, these late/slow shells move in the wake of the external 
shock. Since the circum-burst medium has been swept-up by the external blastwave, the wake should be 
very clean. It means that the movement of the late/slow shells should be nearly constant. So these shells 
must be cold and will not expand laterally. When they finally catch up with the external shock, energy 
injection will only occur at the central region of the GRB remnant, with the involved half opening angle 
(i.e., the energy-injection angle) still being 6'o,n- In Figure 1, we present a schematic illustration of the 
scenario. 



Wide componei 
Narrow componei 




Fig. 1 A schematic illustration of the localized energy-injection process. The outer solid 
curves represent the wide component. The middle solid curves represent the narrow com- 
ponent. Both the wide and the narrow components can expand laterally. The inner dashed 
straight lines represent the trace of the late/slow shells, which coast freely and experience no 
lateral expansion before catching up with the external shock front. The half opening angle of 
the late/slow shells is assumed to equal to that of the narrow jet component, since they are 
most likely be of the same origin. 



3 RAPIDNESS OF THE REBRIGHTENING 

An energy-injection can naturally lead to a rebrightening in the afterglow. In the case of a localized 
energy-injection, the rebrightening should be much steeper as compared to a usual homogeneous energy- 
injection. A few parameters may affect the rapidness of the rebrightening, such as the number density of 
the circum-burst medium the half opening angle of the late/slow shells (6*0,11), the energy fraction of 
electrons with respect to protons (Eo), the energy fraction of magnetic field with respect to protons (cg), 
the power-law index of electron distribution function (p), etc. Analytically, the evolution of the radius 
(i?) and Lorentz factor (7) of the external shock is i? oc n~^/''t^/*, 7 oc ?i^^/^t~^/^ dHuang et al.l 
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Fig. 2 The effects of the parameter 9o,n (a) and n (b) on the rebrightening caused by a local- 
ized energy-injection. In Panel (a), 6*0. n is in units of radian, and in Panel (b), n is in units of 
cm^^. The values of other parameters are set constant in these calculations and are explained 
in the main text of Section 3. 



Il998h . The angular timescale, i.e. the approximate duration of the rebrightening, can then be derived as 

AT w -Rinj(l — cos^o.iO/c ~ ^inj^o n/2c oc u^^^^Oq ^tj^^ , whcrc i?inj is the radius of the external 
shock at the energy-injection moment of tiny It indicates that the effect of Oq ^ on the rapidness of the 
rebrightening is most striking, and the parameter n only has a minor effect on the rebrightening, while 
the effects of other parameters are very insignificant. 

We have numerically investigated the effects of the involved parameters. Our numerical results are 
consistent with the above analysis. In Fig. 2(a), we illustrate the effect of 0o,n on the afterglow light 
curve. In these calculations, we take Co = 0.1, = 0.01, n = 1 cm~'^, p = 1.9, and evaluated 6*0.11 as 0.02, 
0.035, and 0.06 respectively. It is quite clear that a small energy-injection angle can lead to a very rapid 
rebrightening, and a large injection angle can make the rebrightening very smooth. It strongly hints us 
that we could use the rapidness of the observed rebrightening to derive the energy-injection angle. 

As a comparison, we illustrate in Fig. 2(b) the effect of n on the afterglow light curve. Again we 
have taken ec = 0.1, = 0.01, p = 1.9, and 6*0.11 = 0.035, but let n vary between 0.5 cm^"^ and 2 cm^'^. 
We see that the effect of n on the rapidness of the rebrightening is much weaker. 

In short, we find that the rapidness of the rebrightening is primarily determined by the energy- 
injection angle (i.e. the initial beaming angle). We suggest that we could use the observed afterglow 
rebrightenings to measure the beaming angle of GRBs. This is a novel and hopeful method to constrain 
the collimation degree of GRBs, independent to the conventional jet-break timing method. 



4 APPLICATION TO GRB 030329 

In this section, we use the localized energy-injection model to explain the multiband afterglow light 
curves of GRB 030329. Special attention will be paid to th e remarkable r ebrigh tening at t ~ 1.2 x 10^ 



X (12004. However, the original 

data include the contribution from the host galaxy dGorosabel et al.ll2005 ) and the underlying supernova 



s. The observed data of optical afterglow are pro vided bv iLipkin et al 



JZeh et al. 2004, 2005). They are also contaminated by the extinction of our Galaxy (Schlege l et al.l 
fl998h . The pure optical afterglow light curves of GRB 030329 are available only when these contamina- 
tions are corrected for. We will use the se pure afterglow lig ht curves as the final template. The observed 
date of radio afterglow can be found in lBerger et al.l (l2003l) . The effect of synchrotron self-absorption is 
considered in our calculations. 

To fit the overa ll multiband afterglow of GRB 030329, we basically need a two-component jet model 
(iBerger et al.ll2003b . The central narrow jet component can account for the prompt GRB emission and 
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the early afterglow, while the wide jet component can account for the late afterglow. The observed 
multiple rebrightenings are then produced by a few energy-injections. In our framework, the energy- 
injection is caused by the late/slow shells ejected from the central engine. These shells should have the 
same origin as that of those early shells which generate the prompt GRB emission via internal shocks, 
except that they are ejected in a slightly late period and with much lower velocities. It is then natural 
to assume that the initial half opening angle of these shells equal to that of the early shells (i.e., the 
narrow jet component). The collision between a late sh ell and the external shock is very quick. We can 
approximate it as an instant energy-iniection proce s s jRees & Meszaroslll998l: Ik umar & Piranll2000 l : 
ISari & Mesz^l2000t|piran et al.i2004l:lKing et aLlbOOslk 



To fit the observed multiple rebrightenings, we find that at least 5 energy-injections are needed. 
After correctting for the cosmological redshift, the injection time (and the amount of energy needed) is 
1.1 X 10-* s (0.16 £;„,iso), 1-0 X 10^ s (0.15 E^naso), 1-8 X IQS s (0.09 £;„,iso), 2.2 x 10^ s (0.04 £;„,iso), 
and 3.6 x 10^ s (0.06 E'n.iso), respectively, where E'n.iso is the isotropic -equivalent kinetic energy of the 
narrow component. In Figs. 3 and 4, we illustrate our best fit to the multiband afterglow of GRB 030329. 
In these calculations, the initial half opening angle, the initial Lorentz factor and the isotropic -equivalent 
kinetic energy of the wide jet component are taken as 0o,w = 0.2, 70. w = 50, and i?w,iso = 3.4 x 10^^ 
ergs. Correspondingly, the parameters of the narrow jet component are 6'o.n = 0.035, 70. n = 200, and 
-En.iso = 2.0 X 10^^ ergs. The common parameters of the two components are ec = 0.1, eg = 0.01, n = I 
cm~'^, p= 1.9. The luminosity distance is taken as _Dl = 0.8 Gpc. Note that we have assumed a p value 
smaller than 2. Alth ough such a small p is not common for other events, it is still possible in GRBs 
(iDai & Chendl200lh . 

In Fig. 3, we see that the overall R-band afterglow light curve can be well explained. Especially, the 
observed rapid rebrightenings can be satisfactorily accounted for. Taking the most obvious rebrightening 
observed at 1.2 x 10^ s as an example, this rebrightening is s o rapid that all previ ous studies failed 
to present a satisfactory explanation. For example, in the study bv lHuang et al.l (|2006|) . since the energy 
was assumed to be homogeneously injected into the whole external shock, the theoretical rebrightening 
is then much slower than observations, and the relative residual is generally on a level of ~ 20% at the 
rebrightening moment. In our current study, the rapidness problem is satisfactorily resolved. The relative 
residual is only about 5% — 10% during the rebrightening (see the bottom panel of Fig. 4). 

Fig. 4 shows that most of the multiband afterglow light curves can be satisfactorily reproduced. 
The most obvious deviation appears in the 4.86 GHz light curve. Our theoretical flux is lower than the 
observed intensity after 2 x 10^ s. Note that at this late stage, the decay of the observed 4.86 GHz 
flux density becomes very slow and the light curve becomes flatter. It is possible that the contribution 
from the host galaxy begins to emerge from this point on. However, it is also possible that the structure 
of the ejecta in GRB 0303 29 were much more complicated than what we have considered here. In fact, 
der Horst et al. I (120051) have argued that additional component is necessary to account for the excess 
of the observed flux at l ate stages. Finally , radio afterglow is subjected to scintillation due to scattering 



of interstellar medium jFrail et alj|2003l) . This may also explain part of the deviations, especially at 
early stages. Here, since our interest is mainly on the earlier optical rebrightenings, we suspend from 
extending our discussion to more of the above details. 

Many other authors have studied the afterglow of GRB 030329 ( iBerger et alj|2003l: iGranot et all 
l2003l : lHuang et al.l2006tlGao & Dail2010l) . It is interesting to note that the derived paramet ers vary from 
one gr oup to another. For exa mple, the suggested bea ming angle of the na rrow jet is 0.05 bv lHuang et al.l 
(I2OO6I) and lGao & Pail (I2OIOI). 0.07 by ?, and 0.09 bv lBerger et"ai](l2003[). The suggested beaming angle 
of the wide jet is 0.3 by Berger et al. (l2003b . and 0.15 bv lHuang et alT i 20061) . In our current study, we 



take 6*0,11 = 0.035 and Oq.w = 0.2. Although these two parameters are again different from the above 
values, they are not unreasonable, and should be acceptable. Similarly, our other parameters are also in 
reasonable ranges. 
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Fig. 3 Our best fit to the R-band afterglow light curve of GRB 030329, using the local- 
ized energy-iniecti on model. Observed data points correspond to pure afterglow emission 
dLipkin et al.l2004|) . The solid line corresponds to our best simulation with five major energy- 
injection processes which occur at 1.1 x 10^ s, 1.0 x 10^ s, 1.8 x 10^ s, 2.2 x 10^ s, and 
3.6 X 10^ s respectively. The bottom panel gives the relative residual of the fit. See the main 
text in Section 4 for more details of this calculation. 



5 DISCUSSION AND CONCLUSION 

The rapid rebrightening is an interesting feature of the afterglow of GRB 030329. Conventional ho- 
mogeneous energy-injection model cannot account for the rapidness of the observed rebrightenings. In 
this study, it is clearly shown that the localized energy-injection model can satisfactorily resolve the 
rapidness problem. By using a two-component jet model, together with 5 localized energy-injections, 
we can satisfactorily reproduce the overall and multiband afterglow light curves of GRB 030329. In our 
scenario, the narrow component gives birth to the prompt GRB emission and the early afterglow, while 
the wide component contributes significantly to the late afterglow. The energy-injections are induced by 
5 late shells, which should be ejected by the central engine via the same or similar mechanism as the 
narrow jet component, but at a slightly later period and with much lower velocities. These shells coast 
freely in the empty environment after the external shock. They keep to be cold and do not experience 
any lateral expansion until they catch up with the external shock and inject their energies into the GRB 
remnant. The energy-injection is a prompt process, and most importantly, is a localized collision. It can 
lead to a very rapid increase in the afterglow brightness, barely subjected to the smearing of the equal 
arrival time surface effect. 

The initial beaming angle is an important parameter of GRBs . It can provide important clues for the 
central engine dCheng & Lull200ll:IZhangl2007HGao & l5aill20I0h . According to our best fit, the energy- 
injection angle is about 0.035 in the case of GRB 030329. In our modeling, this angle also equals to the 
half opening angle of the narrow jet component. We also showed that the rapidness of the rebrightening 
is mainly controlled by the energy-injection angle. Other parameters, such as e^, eg, p, and n only 
have minor or negligible effect. So, we suggest that the rapidness of the observed rebrightenings can 
basically be used to derive a measure of the beaming angle of GRBs. This is a novel and hopeful method, 
independent to the conventional jet-break timing method. Comparing with the jet-break timing method, 
the advantage of our method is that it does not rely on the assumption of the speed of lateral expansion, 
which itself is quite unclear currently. The impact of another factor, i.e. the density of the circum-burst 
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Fig. 4 Our best fit to the multiband afterglow light curves of GRB 030329 using the localized 
energy-iniection model. Panel (a) co rresponds to optical afterglow. The observed data points 
are derived from lLipkin et alj ( 2004t). P a nel (b ) corresponds to radio afterglow. The observed 
data points are taken from Berger et"ai] ( l2003b . Note that in Panel (a), the I and V band light 
curves have been shifted by —1.2 and +1.2 magnitude respectively, for clarity. See the main 
text in Section 4 for more details of this calculation. 



environment, is also very week in our method. However, the main restriction of our method is that it is 
applicable only when rapid rebrightenings due to energy -injections were observed in the afterglow. 

In our modeling of the five observed rebrightenings of GRB 030329, the beaming angles of the 5 
energy-injections are set to be constant. The fitting results are good. It means that the ejecting angle of 
the central engine might not vary significantly over the active period. However, in the future, it is possible 
that we could find new examples which need to vary the energy-injection angle. Such information will 
also provide interesting constraints on the central engine. In our model, we have assumed that the energy- 
injection angle equals to the beaming angle of the narrow jet that produces the prompt GRB. Although 
this assumption is a natural hypothesis, it is still possible that in reality they might not equal, since the 
real central engine may be much more complicated as compared with our imagination. Even in this case, 
the energy-injection angle derived from observations still can give useful hints on the characteristics of 
the central engine. 

In our current study, in order to get an acceptable fit to the overall and multiband afterglow of GRB 
030329, we have adopted a two-component jet geometry. However, note that such a two-component 
configuration is not always necessary for localized energy-injection. For example, in other GRBs, it is 
quite possible that the prompt GRB jet may only have one component. This component can expand 
laterally during its deceleration, still leaving a relatively clean environment behind it. So, late shells 
will still coast freely after the external shock. They keep to be cold and will not experience sideways 
expansion before the final energy-injection. 

For GRB 030329, in order to get an acceptable fit to the overall and multiband afterglow, we have 
adopted a p value that is less than 2. However, as stated in section 3, p does not influence the rapidness 
of the rebrightening. So p < 2 is not necessary for the localized energy-injection model. In other GRBs, 
different p can be adopted according to observations. 

In many GRBs, multiple optical and/or X-ray flares are observed to be superposed on the early 
afterglow light curves. The rising and falling of these flares are even more rapid than that in the case 
of GRB 030329. Although the nature of these flares is still largely unknown, we suggest that useful 
information on the beaming angle of these flares can also be derived by fitting the rising and falling 
profiles. 
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